High-spin intruder states in the fp shell nuclei and isoscalar proton-neutron 

correlations 
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We perform a systematic shell-model and mean-field study of fully-aligned, high-spin fy^ seniority 
isomers and d~^^f"^^ intruder states in the A~44 nuclei from the lower- /p shell. The shell-model 
calculations are performed in the full sdfp configuration space allowing Ip-lh cross-shell excitations. 
The self-consistent mean-field calculations are based on the Hartree-Fock approach with the Skyrme 
energy density functional that reproduces empirical Landau parameters. While there is a nice 
agreement between experimental and theoretical relative energies of fully-aligned states in N>Z 
nuclei, this is no longer the case for the N—Z systems. The remaining deviation from the data is 
attributed to the isoscalar proton-neutron correlations. It is also demonstrated that the Coulomb 
corrections at high spins noticeably depend on the choice of the energy density functional. 

PACS numbers: 21.10.Pc, 21.10.Sf, 21.60.Cs, 21.60.Jz, 23.20.Lv, 27.40.-fz 



The triumph of the nuclear shell model (SM) in the 
sc?/p-shell region P, Q, 0, Q proves that many spectro- 
scopic properties of those fairly heavy nuclei can be well 
accounted for by an effective two-body G-matrix aug- 
mented by the monopole corrections. In spite of the suc- 
cess of the SM description, there are still many open ques- 
tions and challenges in this region of the nuclear chart 
that offer many opportunities for new physics. In partic- 
ular, studies of mirror-symmetric nuclei and precise mea- 
surements of the Coulomb energy displacement shed light 
on isospin breaking effects 0,|6[. Another frontier is in- 
vestigations of unnatural-parity intruder states in v4~44 
nuclei from the lower- /p shell associated with cross-shell 
excitations across the N=Z=20 magic gap that give rise 
to shape coexistence effects and emergence of collective 
rotational excitations 0, H, 01 ■ 

From a theoretical standpoint, the /p-shell nuclei are 
particularly good candidates to study the competition 
between collective and single-particle excitations. Since 
the associated configuration spaces are not prohibitively 
large for SM calculations, and, at the same time, the 
number of valence particles (and holes) is large enough 
to create substantial collectivity, these systems form a 
crucial playground to confront the spherical SM with col- 
lective approaches based on the mean-field (MF) theory 
|lOl ITU . The diversity of nuclear structure phenomena, 
rich amount of spectroscopic data collected, and a possi- 
bility of direct SM verification of MF calculations in the 
fp region, offer a unique opportunity for fine tuning of 
the underlying nuclear energy density functional (EDF). 

Recently, a systematic MF analysis of maximum-spin 
states (also referred to as terminating states or senior- 
ity isomers) has been performed within the Skyrme- 
Hartree-Fock (SHF) approach [H [3 for the [fy^li^,^ 
and configurations (n denotes the num- 

ber of valence particles outside the ^"Ca core). Those 



fully-aligned states, experimentally known in a number 
of 20 < Z < iV < 24 nuclei, have fairly simple SM con- 
figurations, and they provide an excellent testing ground 
for the time-odd densities and fields that appear in the 
MF description. In this context, the energy difference be- 
tween the excitation energies of the terminating states. 



AE = E{[d-;j^+'] 



J-i?([/7/2]w), (1) 



is a sensitive probe of time-odd spin couplings and the 
strength of the spin-orbit term in the EDF. In partic- 
ular, it was demonstrated [T^ IT^ that by constraining 
the Skyrme EDF to the empirical spin-isospin Landau pa- 
rameters and by slightly reducing the spin-orbit strength, 
good agreement with the data could be obtained. This 
result, based on high-spin data for terminating states, 
is consistent with conclusions of previous works 0, 0| 
based on different theoretical methodology and experi- 
mental input (such as giant resonances, beta decays, and 
moments of inertia) . 

The MF studies of Refs. HI ^ rely on the assump- 
tion that the terminating states are almost ideal exam- 
ples of unperturbed single-particle motion, which further 
implies that AE, unlike the absolute excitation energies 
■^(M3/2/7/2^]/™ax) and £^([/7/2]/,„„J, mainly depends on 
properties of the underlying MF. In particular, AE de- 
pends on the energy of the cross-shell excitation and 
symmetry-breaking effects. The main objective of the 
present work is to (i) study the role of dynamical cor- 
relations on AE, and (ii) investigate the origin of large 
deviations between MF results and experimental data for 
N—Z nuclei. For this purpose, we carry SM calcula- 
tions. Preliminary results of this analysis were published 
in Refs. 18, .19J. 

Our SM calculations were carried out using the code 
ANTOINE in the sdfp configuration space limited 
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FIG. 1: Experimental (gray bars) and shell model (black bars) 
excitation energies of maximum-spin states of /"^j (top) and 
(bottom) configurations in N = Z (left) and A'' > 
Z (right) /p-shell nuclei. Experimental data are taken from 
Refs. 



to Ip-lh cross-shell excitation from the sd shell to the fp 
shell. In the /p-shell SM space we took the FPD6 inter- 
action [201. The remaining matrix elements are those of 
Ref . j2ll|. As compared to the earlier work ) the mass 
scaling of the SM matrix elements was done here consis- 
tently, thus reducing the sd interaction channel by ~4%. 
As seen in Fig. excellent agreement was obtained be- 
tween the SM and experiment for the absolute excitation 
energies of terminating states for both E{[fy^]i^^^) and 

^([^3/2/7/2^]^">»x ) configurations. 

The calculated SM and SHE energy differences AE are 
shown in Fig. |2 relative to experimental values. We note 
that while Fig.nsuggests a similar level of agreement be- 
tween experiment and the SM in N=Z and N>Z nuclei, 
the energy differences tell a different story. Indeed, in 
N>Z nuclei the SM systematically overestimates the ex- 
perimental data by '^280 keV. On the contrary, in N—Z 
nuclei the SM systematically underestimates the data by 
~410keV. This clearly suggests that important corre- 
lations related to isospin and cross-shell excitations are 
missing in the present SM implementation. 

The SM results are further compared to the SHF calcu- 
lations based on the SkO parameterization slightly 
modified along the prescription given in Refs. (12, 17] . 
Without entering into details, we recall that the modifi- 
cations concern coupling constants related to the time- 
odd spin fields and Cf^^s-As where t = 0, 1 labels 
isoscalar and isovector terms, respectively. Moreover, the 
strength of the spin-orbit interaction was reduced by 5% 
compared to the original SkO value. 

In contrast to the SM, the SHF underestimates ex- 
perimental values of AE in N>Z nuclei by '^200 keV 
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FIG. 2: Difference AEth — AE^xp between experimental and 
theoretical values of AE Q in A~44 mass region. Dots de- 
note the SM results. Circles denote the SHF results based on 
the modified SkO parameterization (see text). The SHF cal- 
culations for the [dg^j/T/t^l^moa: intruders in N—Z nuclei yield 
two nearly degenerate states associated with proton (vTph) and 
neutron (vph) cross-shell excitations. As shown in the inset, 
the physical T=0 state in the laboratory frame is shifted down 
in energy by SEt (isospin correlation energy) . Squares denote 
the SHF results for N=Z nuclei with the isospin correction 
added. The SHF results were shifted by 480 keV in order to 
facilitate the comparison with SM. 



giving rise to an average offset of ^^480 keV between the 
two models. In order to facilitate the comparison, this 
average difference was removed by shifting up the HF re- 
sults in Fig. 13 (It is to be noted that an overall shift 
in AE can easily be accounted for by varying the size 
of the N=Z=20 gap in SM or by a changing the mag- 
nitude of the spin-orbit term in SHF.) It is striking to 
see that SHF calculations follow SM results in N>Z nu- 
clei extremely well, reproducing details of isotopic and 
isotonic dependence. This result appears to be fairly 
general. Indeed, as seen in Fig. 13 similar agreement 
was obtained for SHF calculations based on the SLy4 
parameterization modified according to Ref. [T^ . 
These results strongly support our assumption that the 
maximally-aligned states in N>Z nuclei are excellent ex- 
amples of an almost unperturbed single-particle motion 
and that dynamical correlations present in these states 
do not exhibit any distinct particle number dependence. 

The difference in the SHF description of N>Z and 
N—Z nuclei seen in Fig. [3 can be partly explained in 
terms of the spontaneous breaking of isobaric symme- 
try in the (^3/2/7/2^ terminating states in N=Z nuclei. 
In the MF picture, those states are not uniquely de- 
fined. Indeed, by making either neutron (j/) or pro- 
ton (tt) ^3/2—^/7/2 Ip-lh excitation, one arrives at two 



nearly degenerate intrinsic states E{[d^J'^f"j\^^]j^ 



.) 



'^([^3/2/7/Vl/max)' ■^hich manifestly violate isobaric 
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FIG. 3: Difference between SM and SHF values of AE. Two 
Skyrme parametrizations are used: SkO (dots) and SLy4 (cir- 
cles), modified according to Ref. As in Fig. H the SHF 
results were shifted by 480 keV. 



symmetry. Indeed, these MF states are not eigenstates 
of isospin. After isospin projection, the T—0 state be- 
comes lower in energy in the laboratory system, as illus- 
trated in the inset of Fig. |21 Due to physical symmetry- 
breaking caused by the Coulomb interaction, the two in- 
trinsic states are slightly split with the proton Ip-lh ex- 
citation being always slightly lower in energy. The reason 
is that proton excitation from the ^3/2 orbit to a more 
extended /■j/2 orbit slightly increases the mean charge 
radius, thus reducing the Coulomb repulsion. 




FIG. 4: Phenomenological estimates of the isospin energy cor- 
rection, SEt{A), due to the restoration of isobaric symme- 
try internally broken in SHF solutions corresponding to the 
['^3/2/7/2^]^m£i^ terminating states in N—Z nuclei. The values 
of 6Et{A) obtained in SkO and SLy4 models are labeled by 
open and filled dots, respectively. 



In order to make comparison to the data, the corre- 
lation energy 6Et due to isospin symmetry-breaking in 
SHF should be estimated. For the purpose of this work, 
we evaluate SEt using the self-consistent isocranking 
[2^ 1. That is, we compute the energy difference between 
the isobaric analogue states at high spin, i.e., SEt = 
^([^^372/7"/^]^;^. = ±1) ~ E{[d-lJ-+']j_;T^ = 
0) using the SHF approach with Coulomb interaction 

switched off. The energy difference AEj^p corrected 
in this way is marked by squares in Fig. |21 

The calculated isospin corrections are depicted in 
Fig. 0] It is interesting to observe that SEx {A) shows a 
surprisingly strong decrease with increasing A. Accord- 
ing to our analysis, this strong particle- number depen- 
dence can be attributed to the time-odd fields, and the 
calculations indicate that this effect can be reduced by 
decreasing the value of the isovector Landau parameter 
Qq. Whether or not this can be used to further constrain 
the value of Qq remains to be studied (see, however, re- 
cent work [25i. i26i |). Coming back to Figs. |21 and |21 it is 
encouraging to see that after approximate isospin sym- 
metry restoration, one obtains AE'gjjp ' « ^E^m ^^^^ 
in N=Z nuclei. Hence, our comparative study strongly 
suggests that correlations of a similar type are missing in 
N=Z nuclei, both in the SM and the SHF approaches. 

Our SM interaction conserves isospin. Consequently, 
the Coulomb correction to AE, 6Ec should be added 
afterwards. The Coulomb correction (including the as- 
sociated isovector polarization) can be calculated self- 
consistently in SHF. Surprisingly, the many-body re- 
sponse against electrostatic polarization appears to be 
strongly sensitive to the isovector part of the EDF. This 
is visualized in Fig. [S] which shows a difference, SEc, be- 
tween the SHF values of AE calculated without (Ai?j^]^) 
and with (AEhf) the Coulomb term. While SEc is very 
small for SLy4, the values calculated in the SkO variant 
are appreciable, SEc ~ 130 keV. The difference can be 
traced back to the fact that these two parameterizations 
differ strongly in the strength of the isovector part of the 
spin-orbit interaction. While in SLy4 the ratio of the 
isovector {Wi) to the isoscalar (Wo) spin-orbit strengths 
equals to the standard value of Wi/Wq = 1/3, SkO is a 
modern parameterization having Wi/Wq ~ —1.3. The 
resulting change in the radial form factor leads to a large 
Coulomb effect at high spin, an effect that is of the same 
order as the measured Coulomb energy differences in fp- 
shell nuclei |^ . Based on our study, the Coulomb interac- 
tion can give rise to an overall displacement of the order 
of 100 keV that very weakly depends on Z and N. 

In summary, the self-consistent SHF analysis of ter- 
minating states in the v4^44 nuclei agrees nicely with 
SM studies, after correcting the former for the isospin- 
breaking effects in N—Z nuclei. For N>Z nuclei, both 
theories provide a good reproduction of experimental 
data. This validates the assumption of previous stud- 
ies [3 E3 regarding the single-particle character of the 
maximally-aligned states. We believe that the origin 
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FIG. 5: Coulomb correction, SEc, to AE calculated in the 
SkO and SLy4 models by performing SHF calculations with- 
out [AE^Ip) and with (AEhf) Coulomb interaction. 



expected to impact the isoscalar channel associated with 
the sd—>fp pair scattering. The single main obstacle that 
prevents us from carrying out calculations in an extended 
space involving 2p-2h, 3p-3h,..., cross-shell transitions 
is the lack of an appropriate effective interaction. On 
the SHF level, while the extended proton-neutron self- 
consistent SHF formalism has been developed [23|, its 
practical implementation is still in an early stage; see, 
e.g., Refs. and references quoted therein. 

Finally, we have demonstrated that the state- 
dependent Coulomb polarization at high spins noticeably 
depends on the choice of the energy density functional, in 
particular its spin-dependent terms. The resulting uncer- 
tainty in the Coulomb energy shift can be as large as the 
measured Coulomb energy displacement. This is likely 
to result in ambiguities when estimating Coulomb effects 
at high spins. 
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